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The fate of R-parity in the context of the minimal supersymmetric standard model is a central 
issue which has profound implications for particle physics and cosmology. In this article we discuss 
the possibility of testing the mechanism responsible for the stability of the lightest supersymmetric 
particle at the Large Hadron Collider (LHC). The simplest theoretical framework where R-parity 
conservation can be explained dynamically allows for two types of B-L models. In the first scenario 
the new Higgses decay mainly into two right-handed neutrinos giving rise to exotic lepton number 
violating signals together with displaced vertices. In the second model one could have peculiar 
channels with multileptons and/or multiphotons in the final state. In both cases, the local B-L gauge 
symmetry is broken at the TeV scale and the discovery of the new Higgs bosons may be possible at 
the LHC. We investigate in detail the production mechanisms for the Higgs bosons relevant for the 
LHC and the key decays which would shed light on how R-parity is conserved. These results may 
help to understand the link between the cold dark matter of the universe and the missing energy that 
could be observed at the LHC if supersymmetry is realized in nature. 
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I. INTRODUCTION 

The main goal of the Large Hadron Colllider (LHC) is to discover the mechanism responsible for elec- 
troweak symmetry breaking in the context of the standard model (SM) or in a new TeV scale theory. The 
minimal supersymmetric standard model (MSSM) is considered as one of the most appealing contenders 
for this new theory. In this context two important cosmological issues can be solved: the matter-antimatter 
asymmetry can be understood through the electroweak baryogenesis mechanism and the cold dark matter 
of the universe candidate may be the lightest supersymmetric particle (LSP). See Ref. CD for a review on 
phenomenological and cosmological aspects of supersymmetry. 

The fate of R-parity in the context of the MSSM is a central issue which has profound implications for 
particle physics and cosmology. R-parity is defined as R = (—1) 2S M, where S and M = (— 1) 3 ( B ~ L ) 
are the spin and matter parity, respectively. Here B and L stand for Baryon and Lepton number. The 
possible implications of the conservation or violation of this discrete symmetry have been studied quite 
intensively in the last 30 years by many experts. See for example Refs. (2H6|. However, there are only a 
few phenomenological studies of theories which dynamically explain the origin of R-parity. Recently, we 
initiated such a study in Ref. [7] and extend its scope in this article. 

The simplest way to understand the state of R-parity is in the context of a B-L extension of the MSSM, 
where matter parity is just a subgroup of the new abelian symmetry, U(1)b-l- These theories are quite 
simple because only three copies of right-handed neutrinos are needed for an anomaly free theory. Re- 
cently, it was noticed that the minimal B-L model violates R-parity (H, a scenario further motivated by 
string theory [91] . However, since only experiments will reveal the validity of this symmetry, it is important 
to understand the second possibility as well, i.e. the dynamical conservation of R-parity. This is especially 
crucial because observation of missing energy signals at the LHC do not necessarily bare cosmological sig- 
nificance. Therefore, observing both missing energy and the signals discussed in this paper could increase 
the connection of dark matter to missing energy. 

In the simplest framework for dynamical R-parity conservation, B-L is broken at the TeV scale making 
the model testable at the LHC. We discuss the prospects for testing the mechanism for the stability of the 
LSP in two different models which fit in this framework. Our key findings center around the properties of 
the B-L Higgs which can decay into two right-handed neutrinos in the first model and into two sfermions 
in the second case. The final states in the former case are especially interesting since even though R-parity 
is conserved, the final states can violate lepton number. Furthermore the right-handed neutrinos are long- 
lived giving rise to up to four displaced vertices. The main production channels for the Higgses at the LHC 
are investigated in detail and we discuss all possible signals one could use to test the theory of R-parity 
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conservation. 

This work is organized as follows: In Section II we briefly summarize the main implications from 
R-parity conservation or violation. The simplest theoretical frameworks for R-parity conservation are dis- 
cussed in Section III. In Section IV we discuss the decays of the Zbl gauge boson including the effects of 
supersymmetric particles. All production mechanisms at the LHC for the B-L Higgses are investigated in 
Section V. The decays of the physical Higgses are discussed in Section VI, while in Section VII we study 
the most generic signals coming from R-parity conservation. Finally, we summarize our results in Section 
VIII. 

II. SUPERSYMMETRY, R-PARITY AND THE LHC 

The signals indicating a discovery of low scale supersymmetry (SUSY) at the LHC depend on the 
conservation or violation of R-parity. In fact, both the cosmological and phenomenological aspects of the 
MSSM crucially depend on this. It is well-known that one has the following predictions: 

• R-Parity Conservation: SUSY particles are produced in pairs and typically decay via long decay 
chains with multijets, multileptons and missing energy. The latter is due to the LSP, which is stable. 
Detecting missing energy is then a direct evidence for SUSY dark matter. However, while the LSP 
may be stable on collider scales, its stability on cosmological scales is not assured. If the mechanism 
for the LSP stability (R-parity conservation) is also tested it can shed further light on this issue Q. 

• R-Parity Breaking: One can have single production of supersymmetric particles and possible obser- 
vation of lepton and/or baryon violation at the LHC. See Ref. Q for a review and Ref. lfl0l - fT3l for 
recent studies. Lepton number violation stems from non-vanishing couplings of the type LH U , LLe c 
or QLd c , while the presence of u c d c d c lead to baryon number violation. However, the presence of 
both lepton and baryon number violating terms together would lead to catastrophic proton decay 0. 

In general it is easier to discover SUSY at the LHC if R-parity is broken since SUSY particles decay 
to SM final states instead of missing energy, except for the SM neutrinos. In models with spontaneous 
R-parity breaking through the vacuum expectation value of the right-handed sneutrinos ll3l[8l [T4l[T5l . 
only the bilinear term LH U from above exists at the renormalizable level. Furthermore, it is important 
to note that even when R-parity is broken, the gravitino can still be a good dark matter candidate iPToll . 
We postpone discussing the LHC testability of the theories with spontaneously broken R-parity to a 
later article. 
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If SUSY is discovered at the LHC with missing energy, a possible next step is to test the mechanism 
responsible for R-parity conservation. In the simplest case of a gauged B-L symmetry, which we will pursue 
here, the following items should be searched for: 

• The new neutral gauge boson, Zbl, associated with the local B-L symmetry. For a review on Z 
gauge bosons see Ref. ifTTl . See also Ref. |[T8ll . 

• The right handed neutrinos necessary for an anomaly free gauged B-L theory and study their decays. 
One possibility is through the production mechanism, pp — > Z* BL — > NN. See for example Ref. (19\ - 
l2Tj for a detailed study. 

• Identify the properties of the Higgses responsible for breaking B-L. As will be discussed later, these 
have different relationships to the different LSPs (potential dark matter candidates) and so studying 
their properties may also help to identify the dark matter candidate. 

There are several studies on the discovery of the first two points: Z' gauge bosons and right-handed 
neutrinos at the LHC. However, the properties of the SUSY Higgs bosons responsible for the conservation 
of R-parity have not been studied, except in Ref. Q, which expand upon here by studying the Higgs 
production and decay in more detail. 



The simple B-L extension of the MSSM has two different incarnations which carry a mechanism for 
dynamically conserving R-parity. Before addressing these, we briefly review the status of R-parity in the 



III. THEORETICAL FRAMEWORK FOR R-PARITY CONSERVATION 



MSSM. 



As it is well-known, the superpotential of the MSSM is given by 



Wmssm = WrpC + WrpV 



(1) 



where Wr p c is the R-parity conserving part 



Wr p c 



Y u QH u u c + Y d Q H d d c + Y e LH d e c + fi H u H d 



(2) 



and 



W RpV = eL H u + XL L e c + X'Q L d c + X"u c d c d c , 



(3) 
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contains the R-parity violating terms. Gauging B-L forbids the terms in Eq. which all violate B-L by 
one unit. The most straightforward possibility for the new gauge group is 

SU(3)c®SU(2) l ®U(1)y®U(1) B -l (4) 

Since three copies of right-handed neutrinos are needed to cancel linear and cubic B-L anomalies, the most 
general superpotential becomes 

Wb-L = W RV C + Y V LH U V C + Wextra, (5) 

where the last term is model dependent. The particle content and its charge under Eq. ([4]) is that of the 
MSSM: 



Q T = [u, d) ~ (3, 2, 1/6, 1/3), u c ~ (3, 1, -2/3, -1/3), d c ~ (3, 1, 1/3, -1/3), 

L T = (0,e) -(1,2,-1/2,-1), e c ~ (1,1, 1,1), (6) 

^=(3+^)^(1,2,1/2,0), Hj= (H° d ,H-) -(1,2,-1/2,0), 



plus the right-handed neutrinos: 

v c ~ (1,1,0,1). (7) 

The only remaining sector left to specify is the Higgs content which serves to break U(1)b-l an d a l so 
governs the dynamical conservation of R-parity. Here we will introduce two possibilities within the simple 
framework of adding only a vector-like pair of Higgses. In general, we will refer to these Higgses as 

| ~ (1,1,0,^) 0~ (1,1,0,-^). (8) 

• Model I (r/^ = 2): Here we dub the Higgses X,X ~ (1,1,0, ±2). The extra term in the above 
superpotential reads as: 

W<2ra = VXXX + fu c V c X. (9) 



Once the Higgses acquire a VEV, the second term above induces a Majorana mass term for the 
right-handed neutrinos making the neutrinos Majorana fermions. Furthermore, the new Higgses can 
decay at tree level into two right-handed neutrinos. Recently, it was noted that radiative symmetry 
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breaking via the / Yukawa coupling dictates that in the majority of the parameter space R-parity 
is spontaneously broken ll22l . In this paper, we do not subscribe to any high-scale scenario and 
simply assume that the Higgses acquire an R-parity conserving VEV and then study their signals at 
the LHC. Interestingly enough, even though R-parity is conserved, lepton number is still broken and 
could manifest itself in the form of same-sign leptonic final states. For lepton flavor violating rare 
leptonic decays, see ||23l . For the study of other aspects of this model see Ref. l|4l |T9l. 

• Model II (rjcf, = ^f+l)' While it is well-known that Higgs bosons with even B-L charge which 
acquire a VEV conserve R-parity (4), we supplement this by noting that 2pj (2q + 1) with p and q 
integers also conserves R-parity. This includes r\$ = 4, 2/3 and 4/3 for example. This model has 
not been studied before and has distinctly different Higgs physics from Model I. We term the Higgses 
in this case S, S ~ (1,1,0, ±775) and the extra term in the superpotential is simply the mass term: 

Wii2a = VsSS. (10) 

Neutrinos in this case are Dirac fermions and the new physical Higgses do not couple to the MSSM 
superfields at tree level. This scenario is quite interesting because it is so distinct from the previous 
case indicating different signatures for the mechanism responsible for the stability of the LSP and 
give rise to very exotic Higgs signals at the LHC. 

In general models of B — L, such as the Models I and II, kinetic mixing is possible between the Z 
and Zbl- However the mixing is constrained to be less than about 10~ 2 and only plays a role in precision 
physics ll24l . We therefore ignore it for the remainder of this work. 

For the remainder of this section we discuss the details of these two scenarios in a general way. 



A. B-L Symmetry Breaking 

In order to discuss the symmetry breaking in these models in a general way, we use the notation 4>, cf> ~ 
(1, 1, 0, in^). Then, can be X(X) in model I or S(S) in model II. The relevant soft terms for our 
discussions are: 

- Csoft D [a v LH u v c - b+<$ + X -M BL B'B' + h.c. 

+ m\\4>\ 2 + m|H 2 + ml\D c \ 2 + (11) 
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where B' is the B-L gaugino and ... indicates MSSM soft terms. Spontaneous B-L breaking and R-parity 
conservation require the nonzero VEVs for 4> and 4>. Notice that in the above equation one should add the 
trilinear term afD c D c 4> in the case of Model I. Using (<p) = v/y/2 and = v/V^one finds 

t r 1| |2 / 2 i -2\ l - , ^ 2 2,1 2-2 , 2/2 -2\ 2 

^ = gl^l ( w +v ) ~ h <$> vv + 2 m ^ V 2 * V 32 ^ ~ U ' ' ( } 

This form is very similar to that of the MSSM and the derivations that follow mirror those of the MSSM 
with the appropriate replacements. Assuming that the potential is bounded from bellow along the D-flat 
direction leads to the condition: 

2b 4> < 2\ H \ 2 + m\ + v(i\ (13) 

while 

&J> (M 2 + m}) (|/i/ + m|). (14) 
is necessary for a nontrivial minimum. Minimizing with respect v and v one gets 

|^| 2 + rrii - \m\ BL cos 2(3' - b^cot/3' = 0, 

{ (15) 
+ m \+ 2 m2 z BL cos 2/?' -6^ tan /?' = 0, 

with tan fj' = v/v and i r n% BL = Sbl 71 ^ 2 + ^ 2 )/4. These can be recast into the more useful form: 

1 / fn\ tan 2 & — m\ \ 

r*.»=-w-( ' )■ <>« 

6, = ( 2 |„/ + + m?) . (17) 

From here we move on to describe the spectrum details. 

B. Mass Spectrum 

Higgs Bosons: 

The physical Higgs content includes the MSSM Higgses: h, H, A, H^, as well as two extra CP-even 
neutral Higgses, H\ and H2, and one CP-odd Higgs, A^ (X\, X2 and Abl in Model I and Si, S2 and As 
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in Model II). The complex gauge states can be written down in terms of their real components: 

(f) = ^ v + <t>R ^ + ^ I, ^ = 7^ + ^ + (18) 

and related to the physical states through 











ft) 









cos a' sin a' 



sin a' cos a' 





1 






H 









(19) 



(20) 



where G^, is the Goldstone boson associated with breaking B-L and which is eaten by Zbl- The Higgs 
spectrum is completely parameterized by: tan/3', mz BL and 

= sW (21) 
The eigenvalues and the mixing angles in the CP-even neutral Higgs sector read as 



m Hi, a = 2 + T V (m X " m L J 2 + Am z B L m X sin 2 (2/3')) , (22) 



tan 2a' ™L + ™>z BL 



(23) 



tan 2/3' - "J-| Bi ' 

Notice that in the limit, m 2 z BL S> , which will be employed later, the above simplifies to 

m 2 Hl ~ (1 - sin 2 2/3') , (24) 

™ff 2 ~ ™? Zbl sin 2 2/3', (25) 

tan 2/3' ^ 

° ~ _/3 -l + tan 2 2/3' m | M - (26) 

Then, assuming a TeV scale mz BL and small one expects two light Higgses at around the same mass: 
H\ and A&, and a heavy one, H2 close to the Zb l mass. Regardless of the parameter space though, the 
following relationships are observed: mu x < rriA^ and m z BL and m^ 2 > tua^ and mz BL ■ 
Neutrino Sector: 
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The neutrino sector of the two models differs dramatically so we can not discuss the two models generi- 
cally in a worthwhile way. Model II is simple, in this case, neutrinos are Dirac fermions. However, in Model 
I, once the X and X get a VEV a Majorana mass term will be induced for the right-handed neutrinos: 



9BL 



(27) 



noting that / can be diagonalized without loss of generality. This mass in turn triggers the type I seesaw 
mechanism for neutrino masses [25 1: 



m N 



(28) 



As is typical for TeV scale seesaws, Y v ~ 1CT 6 correctly reproduces the neutrino masses. 
Neutralino Sector : 

The neutralino mass matrix, in the basis tp° = ( ^mssm, B', 4>, <j>), reads as: 





M BL 



-MmSSM 



9bl'Jv 



rid, n - 

-9bl -f v g B L -f v 



II® 

2 











(29) 



/ 



where V'sm defines the MSSM neutralinos and .Mmssm the four-by-four MSSM neutralino mass matrix. 
Here, we define the mass eigenstates as in the MSSM 



Xi = 



(30) 



where N diagonalizes the full seven-by-seven neutralino mass matrix and breaks up into block diagonal 
form where the upper four-by-four block diagonalizes the MSSM and the lower three-by-three block diag- 
onalizes the B-L neutralino sector. The eigenstates are labeled with increasing mass so that Xi (x?) i s the 
lightest (heaviest) neutralino, although the lightest B-L neutralino will play a role later so we denote it xbl- 
Sfermion Masses : 

In the sfermion sector, the mass matrices M.\, M. \ and M.\ in the basis f /x, fuj are given by 



m 



Q 



+ K + (5 - 14/) M Z c 2/3 + \D BL 



V2 



Y u fi v d ) 



(o« v u - Y u n v d ) 
I + m 2 u + |M| c 2p s 2 w - \D BL 
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nr, + m\ - (\ - \ s 2 w ) M| c 2/3 + \D BL ^ (Y d fiv u -a d v d ) 

^ (Y d 11 v u - a d v d ) m| c + m 2 d - \ M 2 Z c 2 p - \D BL 
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(31) 



m\ + m 2 e - {\ - sly) M| c 2/ g - D BL ^ (Y e iiv u -a e v d ) 

^ (Y e fiv u -a e v d ) m\ c + m 2 e - M| c 2( g + 

where c 2/ g = cos 2/5, s^/ = sin 6\y and 

Dbl = I qIl^ ( t; 2 - v 2 ) = ^M 2 Zbl cos 2/3'. (32) 

m u , mrf and m e are the respective fermion masses and a u , a d and a e are the trilinear a-terms corresponding 
to the Yukawa couplings Y u , Y d and Y e . Typically, it is assumed that substantial left-right mixing occurs 
only in the third generation. Regardless, the physical states are related to the gauge states by 



cos 6q sin 0q\ ( q 
v Q2 1 \ — sin 6q cos I 




(33) 



where here we are thinking about the squark sector, but, of course, the same thing can be done in the slepton 
sector. 

The left-right mixing in the sneutrino sector is negligible due to the small Dirac Yukawa couplings 
necessary for the type I seesaw mechanism. The left-handed masses are 

m l L = m L + \ M z cos 2/9 - Dbl- (34) 

In Model I the right-handed sneutrino CP-even and CP-odd states are split by trilinear terms involving 
the B-L Higgses. Remembering that the Yukawa matrix, /, can be chosen to be diagonal without loss of 
generality, the masses of the right-handed sneutrinos are given by 

m\ R = m? c + 2f 2 v 2 + V2a fi v + V2f if x x v + D BL , (35) 



m 



Nj. 



m% + 2f 2 v 2 - V2a ft v - V2 f t fi x v + D BL . (36) 



where i runs over all three generations and repeated indices are not summed. The masses for Model II 
can be recovered from the above by setting /j, a/ s — > and in this case the right-handed sneutrino can be 
treated as a single complex scalar field. 

It is important to reemphasize that in this context R-parity conservation and therefore the stability of 
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the LSP is a direct consequence of the breaking of B-L via 4> and 0. The properties of these fields are the 
crucial ingredient for testing this mechanism. These can give rise to unique signals: lepton number violating 
in Model I (despite the conservation of R-parity) and multi-leptons and/or multi-photons in Model II. 



IV. DECAYS OF THE Z BL NEUTRAL GAUGE BOSON 

The discovery of a new B-L gauge boson at the LHC is crucial to establish the existence of a new 
abelian gauge symmetry and to test the mechanism responsible for R-parity conservation or violation in 
the supersymmetric case. In this section we discuss the main features of the Zbl boson decays in order to 
understand the impact of the supersymmetric particles on the total width. Furthermore, Zbl can decay into 
the Higgses Abl and X\, a decay that does not exist in the minimal non-SUSY B-L model (since there is 
no Abl)- As we will discuss later, this decay width also enters into the Higgs pair production cross section, 
pp — > Z* BL — > H\A$, which is an important channel for the discovery of these fields. 

The current bounds on Zbl from LEP II are commonly quoted as mz BL / Qbl > 6 TeV [26] but since 
our covariant derivative is defined as = — \ gsL n<j> Zbl> tne re l evant bound here is 

>3TeV. (37) 
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In what follows, we will simply take this upper limit as an equality. 

The Zbl boson can decay into a pair of fermions, light or heavy neutrinos, sfermions, or into a pair of 
two new Higgs boson. The partial widths for the decay into particles Pi , P2 of masses mi , 771,2 are given by, 



T(Z BL -> P1P2) = rrz ~ \M(Z BL -> P^)? 

167T777z Br , 



(777,1 + fn2) 2 \ j ( m i — fn-i) 2 



1- : " 1 



2 11 2 

ml / \ 777^ 

^BL I \ ^B 



L 



(38) 



where the squared matrix elements follow from the Feynman rules in the Appendix: 



\M{Z BL ^ hfi)\ 2 = \c f (^n f BL ) 2 m 2 ZBL |\ + -J^J , ft = u ,d,c, s,b,t,e, ^r; (39) 
\M{Z BL -+ vm)\ 2 = \ {^n BL ) 2 ml BLl (40) 
\M(Z BL -+ NN)\ 2 = I (cfn B tf m \ BL (l - , (41) 
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M{Z BL ^f a m 



1 (9BL f \ 2 2 

3 C fV~2~ BL ) mZBL 

V 



( 2m 2 - + 2m 2 - 



fa 
2 

Zbl 

fafg = 



+ 



2 2 
771 - — 771 - 

Ja J/3 

4 



in 



(42) 



Zbl 



\M(Z BL -> X^ Bi )| 2 = i (^n* L 



rn 



^ 2m\ +2m 2 A 



2~- 



3 V 2 
x cos 2 (/3' -a') 

- — —Tlr, 



Zbl 



1 



+ 



2 2 



V 



n? 



Zbl 



rn 



(43) 



Zbl 



m 2 + + QmiTUj 



\ 



2m 2 

Zbl 



3 V 2 



Zbl 



2 2 
nt; — 

2m 4 7 



(44) 



Here, % is a generation index, cj are color factors (c qi = 3, q 4 = 1) and ZJf are the unitary sfermion 



mixing matrices introduced in Eq. p3]> and iVjj are the neutralino mixing matrices defined in Eq. (30). 



Using the above expressions we show the branching ratios of Zbl in Fig. [T] In order to simplify our 
analysis, we choose the masses of the three right-handed neutrinos, =95 GeV. We consider one light 
squark, = 150 GeV, and one light slepton, = 150 GeV. All other eleven squarks, five charged and 



six neutral sleptons are heavy (including the three right-handed sneutrino), m„ 



TTlr 



1 TeV All mixing 



angles in the sfermion sector are set to zero for simplicity. The masses of the new neutralinos are determined 
by \ix and Mbl, both are taken here to be 150 GeV. Only the lightest state contributes, while the heavier 
ones have masses very close to mz BL and give negligible or zero contributions. Notice that the numerical 
results are shown for the model I, where the Higgses breaking B-L have = ±2. 

Fig.[T]shows that once the Zbl mass is above 2 TeV, the "susy threshold", the decays into superpartners 
can become important. In the scenario considered in Fig. [TJ for mz BL = 3000 GeV we have the following 
approximate leading branching ratios: 

• £ ~ 24 - 4 % 

• £jj~ 13.6% 

• £ uu, £ NN ~ 12.2% 



SUSY ~ 28. 



14 




FIG. 1: Branching ratios of the Zbl boson for m^ BL — 220 GeV and mx x — 200 GeV. The masses of the three 
right-handed neutrinos are m^. = 95 GeV. We consider one light squark, mj = 150 GeV, and one light slepton, 
m-?! = 150 GeV, all other eleven squarks, five charged and six neutral sleptons are heavy (including the three right- 
handed sneutrinos), rriq = m? = 1 TeV. All mixing angles in the sfermion sector are set to zero. The neutralino 
masses are determined by fix and Mbl> both taken here to be 150 GeV. 
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FIG. 2: The total Zbl decay width as a function of the Zbl mass for a a SUSY spectrum with an Abl mass of 220 
GeV, the lightest stop and lightest stau mass of 150 GeV and all other sfermions at 1 TeV in blue and /ix = Mbl = 
150 GeV determine the neutralino masses. In red, for comparison, is the total width for the non-SUSY case. Note that 
while the decay to X\ and Abl is not a SUSY decay, it does not exist in the minimal non-SUSY B-L model. In both 
cases, all three right-handed neutrinos are assumed to be degenerate with a mass of 95 GeV and rtix x — 200 GeV. 

The total decay width of the L is shown in Fig.|2]assuming all three right-handed neutrinos are degenerate 
with a mass of 95 GeV and the maximum value of qbl consistent with LEP II. In order to further investigate 
the impact of the supersymmetric particles on the total decay width, we compare the total width for a given 
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q a i * 

5 





(a) 






FIG. 3: Parton-level Feynman diagrams for the production of the CP-even Higgs Xi at lowest order, via (a) single 
production (gluon-gluon fusion or qq annihilation), (b) AblX% pair production, (c) associated ZblXi production, 
and (d) Zbl boson fusion. For (d), the diagram with crossed external lines is not shown explicitly. The initial state 
particles q, q' can be any of the light-flavor quarks or antiquarks. 



SUSY spectrum (blue line) with the non-SUSY case (red line). For Zbl masses above the SUSY spectrum, 
the decays into supersymmetric particles contribute significantly and the decay widths can have significant 
variation between the two cases. 

The key difference in the above for Model II is the different value of and the Dirac nature of the 



neutrinos. In this case, the Zbl decay width to Dirac neutrinos is simply given by Eq. (39 1. However, the 
main features of the supersymmetric contribution to the branching rations and total width are similar. 



V. PRODUCTION MECHANISMS OF THE B-L HIGGS BOSONS 



The dominant contributions to B-L Higgs production arise from the single CP-even production via gluon 
fusion and pair production of the CP-even and CP-odd Higgses. Subdominant channels are associated 
X^Zbl production and Zbl boson fusion. The corresponding parton-level Feynman diagrams are shown 
in Fig. [3] 

In the following we focus on the production of the lighter of the B-L Higgs bosons, X\, in Model I. 
The production cross sections for the heavier Higgs boson, X%, follow in complete analogy by replacing the 
corresponding couplings, but they are suppressed by the heavy mass and thus do not play an important role 
for our phenomenological studies. The results for Model II follow by scaling with the corresponding B-L 
factor. 



16 



A. Single Production via Gluon Fusion 



Single production is possible at the one-loop level via gluon-gluon fusion, gg X\, where squarks run 
inside the loop, see Fig.[3](a) and qq — >• X\, at one-loop level mediated by a gluino-squark loop, see the last 
graph in Fig. [3] but this contribution is highly suppressed by the light quark masses and we neglect it. Both 
these channels depend on SUSY interactions of gauge coupling strengths between the Higgs and squarks: 
the D-terms. 

The cross sections for the single production can be given in analogue to Higgs boson production within 
the MSSM |j27ll28l . making sure to include only the relevant diagrams from Fig. [3] The cross section is 
related to the decay width of the scalar and at partonic level it is given by 



TT 



5m 



T Xl ^gg8( 1 



Xi 



in 



where s is the partonic center-of-mass energy squared. The decay width can be written as 



-;/.</ 



9 S 
1287T 3 



£4 



Xi 



with Tq = Ami/m Xi in terms of the kinematic function A(r) = — \ r (1 — r/(r)), and 



(45) 



(46) 



f(r) 



arcsm 



_I flog 1+iG=t 



ITT 



T > 1 

r < 1. 



(47) 



The sum in Eq.(|45|) runs over all twelve squark eigenstates and the couplings g~ 1 = g s * Sn are given in 



the appendix. Note that only the diagonal Xqq couplings enter since the gqq couplings preserve gauge and 
mass eigenstate of the squarks. 

The couplings for the two eigenstates of a given squark q a = qi^ only differ by sign. This allows us to 
rewrite Eq.([43]) as: 



(T 99^X 1 



1024vr 




6 1 



m 



Xi 



(48) 



where now the sum runs over the six squark flavors. From this result one can see that in the case of 
degenerate squark masses, = m ? - 2 , the contributions cancel within each squark flavor, due to the 
opposite B-L charges of the left- and right-handed squarks. 



17 



> 

O 800 





















1 
























- 






\ 












<r> 1 fb 






a 


























































o->5fb 




















10 fb 














<r> 1 fb 
















































































c > 


o.i a 


























































0.01 fb 





























10 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
2 _ .. 6 light squarks_ 



pp->X f 14TeV LHC 



t only 




2000 3000 4000 

«iz BL (GeV) 



50 100 150 200 250 300 350 400 450 500 
m x [GeV] 



FIG. 4: Hadronic cross sections for single X\ production at the LHC. One squark is considered to be light, mf = 150 
GeV and squark mixing is neglected. We use the MSTW 2008 LO pdf [29 1 at a central factorization scale fi = mxj /2. 
In the left panel, plot lines of constant cross section in the ui J \ BL —mz BL with fixed tan f3' = 1.5 in black and in white 
are lines of constant X\ mass. In the right panel, the cross section is shown as a function of mx 1 for vtia bl = 1 TeV 
and niz BL — 1.5 TeV. Here we also explore the possibility of six light quarks (dotted) and the effects of changing the 
lightest stop mass (solid). 



At the hadronic level, the cross section is obtained from the partonic one by the convolution, 



1 d rPP 

dr dL " 



TO 



dr 



(49) 



with t = s/s, s being the hadronic center-of-mass energy squared, and tq = m 2 Xi /s is the production 
threshold. The parton luminosities are given by, 



1C AB 
aL -ab 

dr 



1 + s ab 



dx r 

X 



(50) 



where the parton distribution functions (pdfs) / a /^(x, fi) parameterize the probability of finding a parton a 
inside a hadron A with faction x of the hadron momentum at a factorization scale fi. In Fig. [4] numerical 
predictions for the single Higgs production cross section via gluon fusion are given. The cross section 
depends strongly on the supersymmetric spectrum. 

We consider here the conservative case where only one squark is light (m^ = 150 GeV) and all other 
squarks are heavy and degenerate in mass and thus cancel each others contributions. In the left panel we 
show the curves of constant cross section in the vtla bl — fnz BL plane for fixed tan @' = 1.5 in black. Lines 
of constant X\ mass are shown in white. The plot reflects the sharply peaked nature of the A(r) function 
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from Eq. (46) at mx 1 = Zrn^ where the cross section can rise to about 16 fb but then rapidly decreases 
when lowering the mass of X\ due to the function A(j) and when raiding the mass of X\ due to both A(j) 
and the decreased gluon luminosity. 

In Fig. 4 (right panel) the cross section is given as a function of the Higgs mass mx 1 , for fixed input 
parameters m.A BL = 1 TeV and mz BL = 1.5 TeV and one light squark, the stop (solid lines) and shows 
the effects of increasing the stop mass. The cross section is very sensitive to the stop mass and decreases 
quickly for heavier squark masses. To the left of the peak, the suppression is due to the function A(t) while 
on peak and to the right its due to the decreased gluon luminosity. For illustrative purposes we also consider 
the most optimistic case in which one squark of each flavor is light, i.e. six light squarks (dashed line). As 
one can read from Eq.(|48]), the result simply scales by six, the number of light squarks. In this case the cross 
section reaches 10 2 fb. Unfortunately, this production channel strongly depends on the SUSY spectrum and 
therefore does not allow for general predictions to test the mechanism behind R-parity conservation. An 
interesting property of this channel though is that a light Zbl is not necessary for production. This is 
different for the pair production discussed in the next subsection, which does not depend very strongly on 
the SUSY spectrum. 

B. Higgs Pair Production: pp — > X\A B l 

The Higgs pair production mechanism is the most important channel for our study and part of its interests 
stems from the fact that while it is not a SUSY process and is fairly independent of the SUSY spectrum 
(only via the Zbl width), it does not exist in minimal non-SUSY B-L models. The reason is of course 
familiar to SUSY practitioners, namely that SUSY requires vector like pairs of Higgses since these scalar 
fields have corresponding fermionic fields which contribute to the triangle anomalies. Therefore, a minimal 
non-SUSY theory has only one CP-odd scalar which becomes the longitudinal component of the Zbl while 
in SUSY there are two, one of them physical. To our knowledge and we believe for this reason, this process 
has not yet been discussed in the literature. The production process proceeds via 

q{pi) ?(pa) ->■ *i(pa) a bl {pa)- (51) 

The differential partonic cross sections is given by the spin- and color- averaged squared matrix element, 

2 dps( 2 ) 

da q q^ Xl A BL (s) = \M q q-> Xi A BL (s)\ 2g , (52) 



19 



b 1 







pp A BL Xi, 14 TeV LHC - 


- 




m ZaL = 1 TeV 


tan fl 


= 1.5 


\ \ X 


tan p 


= 2.0 




tan p 


= 5.0 





b 1 



pp — > AgL Xi, 14 TeV LHC 

m z „ L = 1 TeV 

« )ZBL = 2TeV 




4(10 
m A B] 



m x , = 100 GeV 
m x , = 200 GeV 

Wv = 400 GeV 

i i i i i 

200 400 



FIG. 5: j4bl and Xi pair production hadronic cross section for 14 TeV center of mass energy at the LHC as a function 
of the mass m,A BL for fixed values of tan j3' (left) and for fixed values of the Higgs mass mx 1 (right). We use the 
MSTW 2008 LO pdf |29| at a central factorization scale /i = (niA BL +mx 1 )/2. The suppression due to the threshold 
m A BL + nriXi = m z BL is apparent in the solid lines of both plots. The key feature of this production mechanism is 
that it is fairly independent of the SUSY spectrum and that it can be large for a sizable part of the parameter space. 

where dPS^ = di/(8irs) is the two-particle phase-space element. The hadronic cross section follows by 
convolution with the parton luminosities, 



f 1 dL w - 

dcr p p^ Xl A BL (s) = ^2 / dT—^-da q q_> Xl A BL (s), 

q=u,c,d,s T ° 



(53) 



with the threshold being to = (mjd + mA BL ) 2 /s. It is convenient to express the matrix element in terms 
of the usual Mandelstam invariants, 



s = (pi+P2) 2 , i=(pi-p 3 ) 2 , u = ( Pl - P4: ) 2 . 
The squared matrix element can then be written as, 



M q q^ Xl A BL {s)\ 
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tit — m 2 A m? v 

A B L Ai 



(8-m% BL )* + m 2 ZBi T 2 ZBL 



cos 



[ft - of) 



(54) 



(55) 



The numerical cross section results for the pair production of A^l and X\ at the LHC at 14 TeV are 
shown in Fig.[5]as a function of the mass rriA BL - We use the MSTW 2008 LO pdf at a central factorization 
scale n = (rriA BL + mxJ/2. In both plots we use gBL = m z BL /(3TeV). In the left panel we plot 
the pair production cross section versus the mass of Abl for three different values of tan j3' for a fixed 
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FIG. 6: Summary of all production channels for X\ production at the LHC for rriA BL — 500 GeV (right). One squark 
is assumed to be light, mf = 150 GeV (no mixing), the heavy neutrino masses are set to m VR = 95 GeV. We use the 
MSTW 2008 LO pdf at a central factorization scale (half of the final state masses). 

Zbl mass of 1 TeV. With these three numbers the entire Higgs sector is fixed and specifically the mass 
of X\ can be calculated at each point and the larger tan/3' the closer X\ and Abl are in mass, with 
m A BL > fnxx- This specific parameterization has the advantage that the coupling associated with this 
cross section, cos (/3' — a 1 ) is relatively constant over the range of vtla bl shown and so the suppression in 
the cross section for increased in m,A BL is due in small part to the kinematics and in larger part to the Zbl 
threshold (mA. BL + mx x = m z BL ) at which point the Zbl becomes off-shell and the cross section loses its 
resonance enhancement. 

In the right panel, we show the pair production cross section versus rriA BL for three different values of 
mx 1 all for two different values of mz BL (mz BL = 1 TeV solid lines and mz BL = 2 TeV dashed). The 
curves start at the point rriA BL = mx 1 since the CP-even Higgs is at most as heavy as the CP-odd Higgs. 
This plot has the advantage of being in terms of the more physical parameter, mx 1 but then in this case, the 
coupling, cos (/?' — a') changes considerably over the range shown and contributes to the decrease in cross 
section with increasing rriA BL , as does the kinematics. 

In general the cross section for Higgs pair production can be sizable in a large fraction of the param- 
eter space. For example, when the mz BL = 1 TeV, niA BL smaller than 500 GeV and tan/3' > 1.5 the 
cross section reaches several tens of femtobarn. Such results are promising for the prospects of testing the 
mechanism for R-parity conservation. 

In order to complete our analysis we compare in Fig. [6] all of the possible production mechanisms, 
including the associate ZblX\ production and Zbl vector boson fusion. The formulas for the latter two 
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processes are given in Appendix A 2 In most of the considered parameter range, the X\Abl pair production 
dominates and the single production (for mz = 150 GeV) is at a similar order of magnitude. The cross 
section for associate Zbl^i production can be large for light particle but drops off quickly for higher 
masses. The Zbl vector boson fusion, being a 2 — > 3 particle process, is suppressed from the kinematics 
and only reaches the 10~ 2 fb level. In terms of testing the mechanism for R-parity conservation the latter 
two channels play a subleading role and we therefore focus only on the single X\ and X\Abl productions, 
which have the best potential of shedding light on the stability of the LSP. In the reminder of the paper we 
will discuss the subsequent decays of the B-L Higgs bosons and the signals for the B-L Higgs production at 
the LHC. 



VI. HIGGS DECAYS AND LEPTON NUMBER VIOLATING DECAYS 

The decays of the Higgses depend heavily on the spectrum. Here, we will assume that masses are such 
that tree-level two-body decays are open and dominate. The decays can, of course, be very different in the 
two models: 

Model I: 

The two-body decays open to X\ are: 

• Xi NN, 

• x l -> 

Since the coupling of the Higgs to right-handed neutrinos is the defining characteristic of Model I, we 
will assume for the rest of the paper that only the first channel is open and that the SUSY decays are 
not kinematically allowed, namely, mx 1 < 2 mLsp- For heavier right-handed neutrinos, three general 
possibilities exist: decay to one right-handed neutrino and one off-shell right-handed neutrino (rnjsf < 
mjj < 2mw), decays into two off-shell right-handed neutrinos (mx L < m^) or decays similar to Model 
II. Off-shell right-handed neutrinos will manifest as missing energy in final states due to their mixing with 
the left-handed neutrinos. 

The CP-odd scalar, Abl, has the following potential two-body decays: 

• A B l -> NN, 

• A B l ->■ %iXj- 
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• A-bl — > ZblXi, 

where the last one would require a heavy Abl outside the reach of the LHC. The two sfermion channel is 
missing here (compared to X\ decays) since it stems from the D-terms in which Abl does not participate. 
Since Abl is heavier than X\ and we have already assumed that the right-handed neutrino channel is 
opened to X\ it will also be opened to Abl and we proceed by assuming that all others are closed so that 
Abl decays 100% to right-handed neutrinos. 

Therefore, under our assumptions here, the relevant signals to study for Model I are the ones due to the 
decays of right-handed neutrinos, which we will conduct in the next subsection. As we will see, these decays 
could lead to lepton number violating signals due to the Majorana nature of the right-handed neutrinos. 

Model II: 

The difference in the second scenario is the lack of the right-handed neutrino-Higgs coupling, thereby 
only leaving: S\ — > XX an d Si ->■/*/ as possible tree-level two-body decays. If these are not accessible, 
one or more of the following decays will dominate: 

• S\ — > 77 — through a slepton and/or squark loop, 

• Si — > gg — through a squark loop, 

• Si — > Z* BL Z* BL — > some combination of leptons and jets (leptons more likely). 

The CP-odd Higgs now has only two possible decays: As — > x\iXj an d As — > ZblSi- Both decays are 
likely to be outside the kinematic range of the light As we are considering so that one of the final state 
particles in each of these will have to be off shell. Since the latter is independent of the SUSY spectrum, we 
will assume it is the dominate decay with the Zbl off-shell so that the following decays are possible 

• A s -> SiZ* BL -> Si^F, 

• A s -> SiZ* BL -»• Sivu. 

• A s -> SiZ* BL -> Sijj. 

Therefore, we will assume these three body decays for As and two body decays for Si, specifically the 
scenarios where the lightest Higgs, Si, decays into two sleptons: 

• Si -> ee*, 

• Si -> v c v c *. 
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The final states will depend on the identity of the LSP (potentially the dark matter candidate of the universe). 
We will consider the following possibilities: neutralino, gravitino and right-handed sneutrino and their 
associated signals. 



A. Heavy Neutrinos Decays 

Higgses decaying mainly into two right-handed neutrinos allow for lepton number violating signals due 
to the subsequent decay of the heavy Majorana right-handed neutrinos. The leading decay channels for the 
three heavy neutrinos, N a , include: 

N a ^efwT, N a ^u £ Z, N a ^u e h k , N a ^l±H*. (56) 



The amplitude for the two first channels are proportional to the mixing between the leptons and heavy 
neutrinos, while the last one is proportional to the Dirac-like Yukawa terms. While decays to all the MSSM 
Higgses are possible, typically, only the lightest MSSM Higgs, h, is light enough for the scenario we 
consider here (mjv a < 500 GeV) and so we will only take it into account. The partial decay widths of the 
heavy Majorana neutrinos iVj are then given by ll2Tll3"0l 

s W^-^WL^-J)". (57) 

F M T = T{Na £±w ^ = jL\ Vla \2 mNa (l-™^) , (58) 

6Z7T y m Na J 

T^Zl = T{Na ViZl) = _£__ m 2 m Z Na ^ _ J™!. } ^ (59) 



W \ m N a 

.2 / _2 N 2 



r » e z T = T{Na ^ ViZt) = _JL_\ Vea f mNa | 1 _ | ^ (6Q) 



3z,ttc w \ m Na 



2 

L W \ " l N a 



Here the leptonic mixing between the SM charged leptons (£ = e, fi, r) and heavy neutrinos (N = 1, 2, 3) 
reads as [21 J: 

V tN = Vpmns my 2 n M N 1/2 , (62) 
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a product of matrices related to the neutrino sector. Vp m ns is the PMNS active neutrino mixing matrix. 
Under the assumption that it is real: 



VpMNS 



C12 C13 C13 S12 S13 

-C23 Sl2 - C12 S13 S 2 3 C i2 C 2 3 ~ Sl2 «13 «23 C13 S 2 3 
\ S12 S23 - Cl2 C 2 3 S13 -C12 S 2 3 - C 2 3 Su S13 C13 C23 / 



(63) 



where Cy = cos 0y and Sjj = sinfljj with < 0^ < tt/2. For our numerical predictions we assume the 
tri-bimaximal ansatz: 



s 12 



1 



'13 



0, 



s 23 



1 



(64) 



The physical neutrino masses are contained ). As it is well-known, there are 

two possible neuuino spectra: 



Normal Hierarchy (NH): m Ul , m V2 
Inverted Hierarchy (IH): m Vl - 



m 2 ^ + Am 2 



21' m ^3 



m\ + |Ara|,J; 



m 2 , + lAm^ 



to 2 , + Am, 2 ,, , rrij, ; 



where [31 



(65) 



7.27 x lCT 5 eV 2 < Am^ < 8.03 x 1(T 5 eV 2 , 
2.17 x 1CT 3 eV 2 < lAmfil < 2.54 x 1CT 3 eV 2 , 



(66) 
(67) 



are the solar and atmospheric mass squared differences, respectively. In this paper, we will only use the 
central value for these masses. Finally, [32] is a complex orthogonal matrix, which conveniently parame- 
terizes the leftover unknown degrees of freedom of the neutrino sector. We shall proceed by assuming O to 
be real. In this case it can be parameterized by three values: 
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(68) 



w l 2 y 



Since the penultimate final states of interest are composed of right-handed neutrinos, investigating their 
decay properties are worthwhile, especially since they depend on the very small neutrino parameter Vi a . In 
Fig.|7]we do this by plotting the decay length in millimeters versus the mass of the right-handed neutrino (N± 
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iVi N 2 N 3 


BR(iVj -»■ e" W+) 
BR(iV 4 //- 
BR(iVi -»• r" W+) 
£BR(A^i/Z) 


31.9% 15.9% 0% 
8.0% 15.9% 23.9% 
8.0% 15.9% 23.9% 
4.3% 4.3% 4.3% 



TABLE I: Branching ratios for the right-handed neutrinos in the special case of tri-bimaximal mixing, central values 
for the squared mass differences and 51 = 1 for degenerate right-handed neutrinos masses of 95 GeV. 

- red, N2 - blue and N3 - black) in the NH (IH) on the left (right). We scan over the unknown parameters: 
lightest left-handed neutrino mass between 10~ 4 eV and 0.4 eV (where the latter is the upper bound from 
cosmology) and < < 1 for i, j = 1..3. The decay length always increases with decreasing lightest 
neutrino mass for all other parameters constant. 

The noteworthy result from Fig. |7]is that for this range of right-handed neutrino masses — a mass range 
chosen to be consistent with a B-L Higgs being produced at the LHC and decay into two on-shell right- 
handed neutrinos — the right-handed neutrinos are long-lived (order of millimeter or above) and their decays 
would exhibit displaced vertices. This is a robust prediction that would lead to spectacular signals and could 
play a major role in distinguishing these channels. We will expand on this in the next section. 

As can be appreciated from the above, the right-handed neutrino decays can be quite different in a 
given neutrino mass spectrum. To simplify our analysis, we will assume the following: = 1 and that 
the right-handed neutrinos are degenerate in mass. This is in addition to our earlier stated assumptions of 
tri-bimaximal mixing and central values for the squared mass differences. Ref. [21] studies the effects of 
varying these values on the decays of the right-handed neutrinos. 

Under these assumptions the branching ratios are straightforward and independent of the mass of the 
lightest neutrino and the neutrino mass hierarchy and are displayed in Table [I] for degenerate right-handed 
neutrino masses of 95 GeV. These branching ratios would change as the right-handed neutrino mass in- 
creases due to the increase in strength of the Z channels and eventually the Higgs channel — kinematically 
not allowed for these masses — and would eventually level off. Clearly, the branching ratios mirror the 
tri-bimaximal mixing due to O = 1. 

These considerations will impact the final states and therefore the signal. We elaborate on this using the 
simplifying neutrino sector assumptions mentioned above and focusing on the single X\ production and the 
pair production of X\ and Abl- In both cases, the Higgses decay into two right-handed neutrinos which 
subsequently decay into two leptons and two heavy vector bosons. The latter will further decay into jets or 
leptons. The signals we will focus on in the next section are the one associated with lepton number violation: 
the two right-handed neutrino decaying into like-sign leptons (muons or electrons) and W bosons, which 
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FIG. 7: The decay length for the right-handed neutrinos (A^ - red, N2 - blue and N3 - black) versus their mass in 
the normal hierarchy (inverted hierarchy) on the left (right). We scan over the following parameters: lightest neutrino 
mass between 10 -4 eV and 0.4 eV and < < 1 for i,j = 1..3. Due to the small right-handed left-handed 
neutrino mixing which facilitate these decays, the right-handed neutrinos can be quite long-lived which would lead to 
displaced vertices. 
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subsequently decay purely into jets. 

VII. SIGNALS AT THE LARGE HADRON COLLIDER 

We are now ready to study the possible signals of dynamical R-parity conservation at the LHC. For 
each model, we will outline the final state signals for both the single and pair productions and conduct our 
analysis on the cross section times branching ratio level only. We will comment on the relevance of the 
background but of course our comments here would be superseded by a more detailed study of these events. 

Model I : 

As a reminder, we will proceed under the assumption that both the CP-even and CP-odd Higgses decay 
only into two right-handed neutrinos. This allows for final states consisting of like-sign leptons, an indicator 
of the Majorana nature of the right-handed neutrinos, as long as the W bosons decay hadronically. In the 
case were these decays are not possible, the Higgses could decay through off-shell right-handed neutrinos 
or in a fashion similar to Higgses in Model II depending on the spectrum. 

Single Production 

pp -> Xi -»• N N -»• efW T efW T -»• efefAj. (69) 

To get a quick naive estimate for the number of events in this channel we do a back of the envelope calcu- 
lation using 

N 2e 4j ~ <J{PP -> Xi) x BR(Xi -»• TViiVi) x 2BR(iVi -> e + W'f x BR(VF -> jjf x C. (70) 

Assuming a large luminosity, C = 100 jb~ x , and a large cross section of 10 fb one obtains naively: 

N 2e 4j ~ 10/6 x (1/3) x 2(3/10) 2 x (6/9) 2 x 100/fT 1 = 27. (71) 

Then, indicating that a significant number of events can occur. The exact number of events, N enem 4j, can 
be calculated taking into account the contributions of all the right-handed neutrinos using the following 
expression: 

Ne n e m Aj = <J{PP -> X X ) X Y, BR (^1 -»• N * N *) X N le m X BR(H^ -> j jf X C, (72 ) 

a=1...3 
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Final State 


Combinatorics 


Signal 


Background 


2 e ± 4 j 


0.038 


62 


6 


e ± fi^ 4 j 


0.030 


50 


12 


2/x±4j 


0.027 


43 


6 



TABLE II: Number of events for the three possible two same-sign leptonic final states (with e and //) for an integrated 
luminosity of 100 fb _1 and a single production cross section of 16.3 fb corresponding to benchmark I with degenerate 
95 GeV right-handed neutrinos. We also display the combinatoric factor associated with the branching ratio of the 
Higgs to right-handed neutrinos, right-handed neutrinos to the specific leptonic final state and W bosons to jets. This 
factor is independent of the cross section or integrated luminosity and multiplies these two numbers when calculating 
number of events. 

where A/" e a n e m is the combinatorical factor for two right-handed neutrinos decaying into e^e^, 

Nl em = 2BR(iV a -> e+W~) x BR(iV a -> e+W~) x — 2 — . (73) 

We choose a benchmark scenario in order to produce more concrete numbers: 
• Benchmark I: 

_ mABL = lTeV, m Xl = 300 GeV, mz BL = 1.5 TeV 

- m,£ =150 GeV, all other sfermions at 1 TeV 

- rriN i = 95 GeV for i = 1..3 

- In this case Up^Xx = 16.3 fb. 

Using these values we display the predicted number of events for the benchmark I in Table II as well as the 
combinatorical factor associated with the branching ratios of the right-handed neutrinos to charged lepton 
final states and W to jets. This number is independent of the cross section and integrated luminosity and 
multiplies both to find the number of events. 

Meanwhile, the SM background to this sort of signal has been studied before in Ref. [21], and was 
found to be dominated by tiW^ production, with a cross section of 4 fb. Using BR(t — > j&fv) ~ 10%, 
BR(l^ =t — > iff) ~ 10% and BR(i — > jjj) ~ 67% we also include an estimate of the number of 
background events in Table II. Two important points are worth noting about the background. The first is 
that the SM background contains missing energy due to the neutrinos Ell . The second and more important 
is that, as we saw in the last section, the right-handed neutrinos travel a distance of order millimeters 
before decaying thereby producing displaced vertices, a further powerful handle on the signal over the 
background. Therefore, using the information about the two displaced vertices in this case one can suppress 
the SM background. In order to understand the reconstruction of these channels one can use the fact that 
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the invariant mass of two jets should be equal to My/, and the invariant mass of two jets and one lepton 
corresponds to the mass of the right-handed neutrinos EH. 

There is also a possible non-SM background from Zbl — > NN, also studied in ET1 . which would of 
course have the same signal. For masses similar to those in benchmark I, the Z' channel will dominate. 
Assuming that the Zbl mass is known from the electron or muon channel, the reconstructed mass of the 
intermediate particle can be used as a handle to differentiate these two channels. For Zbl masses too heavy 
for the LHC, the single Higgs production might dominate and act as a complimentary discovery channel for 
this model. 

Pair Production 

A very important channel is the pair Higgs production through the B — L gauge boson 

pp -> Z BL -> XiAbl NNNN -> efefe^efSj, (74) 

because it does not depend directly on supersymmetric particles masses and allows for a more reliable signal 
for this mechanism stabilizing the LSP without depending on the SUSY spectrum. 

We again perform a naive estimate to understand the predictions for the number of events with four same 
sign leptons and eight jets signal using a cross section of 100 fb and an integrated luminosity of 100 fb : 

N Ae8j ~ <?{PP -> X^bl) x BR(X 1 -> iViJVi) x BR(A BL -> AWi) x 
2BR(iVi e + W-f x BR(W jjf x C 
= 100/6 x (1/3) x (1/3) x 2(3/10) 4 x (6/9) 4 x 100/6" 1 w 4. (75) 

Here we pick a second benchmark scenario: 
• Benchmark Scenario II: 

_ mABL = 220 GeV, m Xl = 200 GeV, m Zl3L = 1 TeV, 

- m£ =150 GeV, mfj = 150 GeV and all other sfermion at 1 TeV 
_ hbl = 150 GeV, M BL = 150 GeV 

- rn/Vj = 95 GeV, for i = 1..3. 

- The cross is a P p^x!A BL = 65.7 fb. 

We display the predicted number of events in Table III for the five possible final states with e and fi leptons. 
We also show the combinatorics factor which takes into account the branching ratios of the Higgses into 
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rinal Mate 


Combinatorics 


Number of Events 


4e ± 8j 


0.00072 


4.8 


3e ± 


0.0012 


7.6 


2e ± 2^8^ 


0.0015 


9.7 


e ± 311* 8 j 


0.00081 


5.3 


4/i±8j 


0.00035 


2.3 



TABLE III: Number of events for the five possible four same-sign leptonic final states for an integrated luminosity 
of lOOftr 1 and a pair production cross section of 65.7 fb corresponding to benchmark II (degenerate 95 GeV right- 
handed neutrinos). We also display the combinatorics factor which combines the branching ratios for the Higgses into 
right-handed neutrinos, right-handed neutrinos to specific leptonic final states and W bosons into jets. This factor 
is independent of the cross section or integrated luminosity and simply multiplies any cross section and integrated 
luminosity to give the total number of events. 

right-handed neutrinos, right-handed neutrinos into leptons and W bosons to jets. This number can be 
multiplied by the cross section and integrated luminosity to yield the number of events. Note that this 
second benchmark is in some ways complimentary to benchmark I with respect to the mass of Abl- 

In this case the main SM background is ttW^ttW^ which has a negligible cross section. It is important 
to mention that in this case one has four displaced vertices making the signal quite special. This does not 
change the fact though that the reconstruction in this case is quite challenging due to the presence of eight 
jets in the final state. Imposing the condition that the invariant mass of two jets, \M(jj) — M\y\ < 15 
GeV EO, can improve the reconstruction process as well as the order millimeter displaced vertices due the 
long lifetimes of the right-handed neutrino. A more detailed study will be considered in a future publication. 

Model II : 

In this section we assume that = 4 and that the CP-even Higgs Si, decays into two sfermions 
while the CP-odd Higgs, As decays into Si and opposite-sign lepton pairs from an off-shell Z* BL (jets 
and neutrinos are also possible). Regardless of these concrete assumption, the signals still depends on the 
SUSY spectrum. We will therefore only highlight some interesting scenarios and finish by addressing some 
alternatives to the two-body sfermion decays. 

• Si -)• ee*: xi astheLSP: 

Here we assume that < mg < ms 1 /2 and that the lightest neutralino is the LSP (not 

necessarily a B-L neutralino) and that e is the NLSP 

Single Production: In the case of the single production one has 



pp -> Si ->■ e* e -)■ e e T xiXi, 



(76) 
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and it yields opposite-sign lepton and missing energy. For benchmark I, where the single production 
cross section is independent of so that the single production cross section is unchanged {(J pp ^Si = 
16.3 fb), and the number of events assuming 100% branching ratios and C = 1 fb -1 is 

N e±eTE ^ss = a pp ^ Sl x BR(Si -)• e* e) x BR(e exi) 2 x C ~ 16. (77) 

Notice that this estimation is naive because one has to assume a large branching ratio for the decays 
into selectrons. The main SM background is the WW, ZZ, and ti production, but cutting on large 
missing EJj} 1 ss one could reduce this background. See Refs. Il33ll341 for recent studies and examples 
of different techniques. 

Pair Production : As we discussed above the pair production can give us a better idea of the cross sec- 
tion without assuming a particular supersymmetric spectrum. In this case one can have the following 
signals: 

pp StAs -> e* e Siefer -> e ± e T e ± e T efe~xiXiXiXi- (78) 

Then, in this case one has three pairs of leptons and missing EJp tss . This cross section does depend 
on the value of n^, and for benchmark II, <r pp ^SiA s = 160 fb when = 4 (note that the Zbl 
width also changes with so that the cross section doesn't simply scale with this parameter). The 
number of events for C = 1 fb -1 can be estimated as 

N 3(e±eT) E ™°° = ^pp^SxAs x BR(Z BL -> e+e~) x £ ~ 40, (79) 

where BR(Zbl — > e t e 7) ~ ^5% m this case. As in the single production case, demanding a large 
missing one should be able to reduce the background which is much less severe since it involves 
more gauge fields or three pairs of top quarks. 

• Gravitino LSP and xi NLSP: 

For which we assume the hierarchy: < < mg < ms 1 /2 and that neutralino decays 

within the detector: xi jG. 

Single Production: 

pp -+ Si -> e* e -> e ± e T -f-fGG, (80) 
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is marked by a pair of opposite-sign leptons, photons and missing E?- The number of events is 
the same as in the previous scenario since we assuming that all branching ratios are 100%. For 
benchmark I and C = 1 fb _1 : 

N e ± e ^ llE rniss = (Tpp-^Sl X £ ~ 16. (81) 

Again, this is a naive estimation of the number of events. It is easy to see that one can satisfy the new 
bounds from CMS on gauge mediation [35]. 

Pair Production : One can have channels with multileptons and multiphotons if one uses the Higgs 
pair production: 



pp -> SxAs ->■ e* e Sieje' -»■ e + e+ ej ej TrryGGGG (82) 
produces multileptons and multiphotons and missing Et- The number of events is given by 

N 3(e±^)4r,E™*° = °pp-+SiA BL x BR(Z BL -> e+e") x C ~ 40, (83) 

for benchmark II, = 4 and 1 fb _1 of integrated luminosity. One can see that these results are in 
agreement with the bounds from CMS |[35l . 

v c as the LSP: 

In this case the right-handed sneutrino can be a dark matter candidate and in principle the Higgs 
S\ can decay mainly into dark matter, while As decays into two leptons and dark matter, As — > 
S\Z* BL — > S\efeJ — > {u c )*u c ef e~. The number of events can be estimated as in the previous 
cases. However, since in order to study these channels one needs to make use of the initial state 
radiation (ISR), we postpone this study for a future publication. For a study on sneutrino dark matter 
in this context see Ref. [361. 



LSP heavier than S\ : 

Here again there is a lot of variability depending on the specific spectrum. We simply refer the reader 



to Table IV for the possible final states in this case and leave the calculations for the number of events 
to a future paper. It is important to mention again that in the case of Model II it is not possible to 
make well-defined predictions for the signals because we do not know the SUSY spectrum. If one 
sticks to a particular SUSY breaking scenario one could see which of these signals are the relevant 
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TABLE IV: Channels with multileptons and multiphotons in Model II when m LSP > m Sl . 



Si Decay 


Single production final state 


Pair production final state 


Si ->• 77 
Si ->• jj 
Si -> ^*f£ 


77 

jj 


^^7777 

'i'ijjjj 

f± pT p± p^F p± pT p± p^F p± pT 
*-i c i t j c j z k c k z l c l c m c m 



ones. In this paper we pointed out the different possibilities and a detailed study is beyond the scope 
of this article. 

VIII. SUMMARY 

The possibility to test the mechanism responsible for the stability of the lightest supersymmetric particle 
at the Large Hadron Collider has been investigated in detail. As it has been discussed in this article, the 
simplest theoretical frameworks where R-parity conservation can be explained dynamically are based on 
B-L gauge symmetry. We discuss two different models and find the following interesting results: 

• In the simplest theoretical frameworks where one can explain dynamically the conservation of Im- 
parity one must have new Higgs bosons which decay mainly into two right-handed neutrinos or into 
two sfermions. 

• We have investigated the production mechanisms and decays of the B-L Higgses. We have found 
that the Higgs pair production mechanism is quite relevant for the testability of the mechanism for 
R-parity conservation, because its predictions are independent of the supersymmetric spectrum. 

• In Model I, where the B-L Higgs couples at tree level to the right-handed neutrinos, one can have 
lepton number violating signals with multileptons and multijets. In this case, if the masses of the new 
Higgses are below 500 GeV one obtains multiple displaced vertices due to the presence of long-lived 
right-handed neutrinos. 

• A new class of models for the dynamical conservation of R-parity has been discussed. In this case 
the new physical Higgses couple only to the sfermions at tree level and the neutrinos are Dirac 
fermions. One finds different exotic signals. However, those channels depend on the supersymmetric 
spectrum. In a simple scenario, such as gauge mediation, one can have channels with multileptons 
and multiphotons. 
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The testability of the mechanism for R-parity conservations may help us to understand the link between 
missing energy at the LHC and the cold dark matter of the universe. 
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Appendix A: Feynman Rules and Cross Sections 
1. Feynman Rules 

• Higgs-squark-squark: 

(Al) 
(A2) 
(A3) 
(A4) 
(A5) 
(A6) 

for any squark q. 
Higgs-fermion-fermion: 

X^iNr- i 2V2 fi cos a' = -2i g BL °° S " mNi , (A7) 

sin p' m Zl3L 

X 2 N i N l : -i2V2f iS ma' = 2ig BL Sma a/ mNi , (A8) 

sin p' mz BL 

A BL NiN t : 2V2 f t cos P' l5 = -2 gEL-*— (A9) 

tan P' mz BL 



Xicftq! : 


■ x 1 


I 

" 6 


9BLmz BL s'm(a + P') cos 2^, 


Xiq 2 q 2 : 






g 9bl mz BL sin(a + /?') cos 20 q, 


Xiqlq 2 ■ 






g 9bl mz BL s'm(a + P') sin 20q, 


X 2 qlqi ■ 


<\ 




g 9bl mz BL cos(a + /?') cos 20 q, 


X 2 q 2 q 2 ■ 


■ x 2 

l 9q 2 92 


i 

~ 6 


9bl m ZBL cos(a' + 0) cos 20 q, 


X 2 qlq 2 : 


ig q\qi 


i 

~ 6 


9BLmz BL cos(a' + /3') sin 2% 
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quark-squark-gluino: 

9 a ^q*a(p)Uk) ■■ ~ i9s(P + kYKp, (A10) 

g a »g bv q* a qp: ig\ (\ a X b + X b \ a ) <T, (All) 

\ /a/3 

g c^~a . _ g3 f abcl ^ (A12) 

fc9iV a: ^3\/2A^(cos^P L + sin^P R ), (A13) 

^92/35° : igayfiXlf} (- sin %P L + cos 0,-P fl ) , (A14) 

where A a are the generators of SU (3), a and (3 represent color and f a bc are the structure constants 
for SU(3). 

ZblZbl4>i- 

z bl z bl<I>i ~ igBLriBLmz BL sm{P' - OL)g lxv , (A15) 

Z BL Z BL<fa '■ ~ igBLnBLmz BL cos(/3' - ot)g lxv , (A16) 



Zbl&Abl'- 



Z^ L <piA BL : gBL 2 cos ( a ' ~ 

(A17) 

Zbl&Abl ■ 9bl~y~ sin(a' - - P2) M > (A18) 



ZblIJ- Here f = u,d, e. 



f 

Z%Jf: -ig B L^^, (A19) 



Zb L ^z/ and Z B lNN: Here N = vr + {vr) c , 



: *^7m75, (A20) 



^7^75, (A21) 
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ZblIIIj: Here f t = df, ef, uf, where i,j = 1, 2 and a = 1, 2, 3. 



JJ£ = " (p! - p 2 ) M (/.^r/, + C^t^) . (A22) 



^BLXiXj 



Z^XiXr- -^bl^In.^.-N^N^ 7m | (! + «%)• (A23) 
2. Cross Sections 

The process pp — >• Z* BL — > X\Zbl is described as associated production or Higgs strahlung. The 
differential partonic cross sections is given by 

2 dps( 2 ) 

da q q^ Xl z BL {s) = \M q q^ Xl z BL (s)\ — ^— , (A24) 
in terms of the matrix element, 

\M q ^ XlZBL (*) | - ^ (^^^ J (^-m| s J2 + m | fli r| si Sm (/? " " } • (A25) 



The hadronic cross section follows by convolution in analogy to Eq.(p3|> with the production threshold being 

to = (raxi +rnz BL ) 2 /s. 

The result for the Zbl boson fusion, q{pi)q'{p2) — > q(P3)Xi(jp4)q' (ps), arises from the diagram shown 
in Fig. [3^d) and the one with crossed external quark lines. In terms of extended Mandelstams, in = 
(pi — Pi) 2 and U2i = (p2 — Pi) 2 , we can write for the squared matrix element, 

\M gg ^ Xiq As)\ 2 = I (^f^) 2 -l flL sm 2 {ff ~ «') 
x i [s 2 + s(iu + u 2 4 - U23) - (ii3 + i M )u23 + m 2 x (n 2 3 - «)] 7^— — s — — s — ^ 

I (^5-^| BL ) 2 (tl3-m| s J2 ^ 

+ [s 2 + s(n 2 4 + *14 - *13) - C"23 + ^24)il3 + m Xl (*13 - s)] 



(*15-™| si ) 2 (M23-"l| BL ) 2 



+ 



3 

-s(s + iu + u 2 a - m Xl ) 



1 



(u 25 - m 2 ZBL )(ii 5 ~ m 2 ZBL )(i l3 - m| s J(n 2 3 - m| B J 



The matrix element does not depend on the electric charge or the flavor of the quarks and at the hadronic 
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level we can just sum over all possible initial states by adding the respective parton densities: 



dapp^ Xi qq'{s) = 

q,q'=u,c,d,s 




M qq ^ Xiqq ,{s)\ 2 ^dVS^\ (A27) 



where dPS^ is the 3-particle phase-space element and the production threshold is tq = m 2 x /s. 
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